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^ ! Abstract 
00 . 

' We present a new technique - photon folding - for imaging in non-focusing tele- 

I scopes. Motivated by the epoch-folding method in timing analysis, the photon fold- 

^ ' ing technique directly provides the statistical significance of signals in images, using 

I projection matrices. The technique is very robust against common imaging problems 

' such as aspect errors and non-uniform backgroimd. The technique provides a deter- 

"q^I ministic recursive algorithm for improving images, which can be implemented on- 

I ' line. Higher-order photon folding technique allows a systematic correction for cod- 

O I ing noises, which is suitable for studying weak sources in the presence of highly vari- 

> ' able strong sources. The technique can be applied to various types of non-focusing 

I telescopes such as coded-aperture optics, rotational collimators or Fourier grid tele- 

' scopes. 
. ^ I kejrwords: coded-aperture system, modulation collimator, Fourier grid telescope, 

^ ' epoch-folding method 



1 Introduction 



Localization or image reconstruction of astrophysical sources at hard X-ray /y-ray en- 
ergies by non-focusing telescopes, such as a coded-aperture system or modulation col- 
limator, has always been a challenging topic. This is largely due to an intrinsically 
low signal-to-noise ratio (SNR) in non-focusing systems compared to focusing instru- 
ments [|lp. Consequently, the next generation hard X-ray telescopes, such as HEFT]^ or 

*High Energy Focusing Telescope (HEFT) is a balloon-borne experiment that will use depth-graded mul- 
tilayer optics and Cadmium Zinc Telluride pixel detectors to image astrophysical sources in the hard X-ray 
(20 - 100 keV) band. 
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Constellation-X/HXl^ push focusing optics up to energies as high as 100 keV [Q, |]. In 
order to study astrophysical sources at energies above ~ 100 keV, however, we still have 
to rely on some sort of non-focusing schemes for the time being 

For example, gamma-ray bursts, one of the long-standing mysteries in astronomy, 
usually have about 90% of their measured flux in the energy range from ~ 100 keV to 
a few MeV. Their secret is expected to be solved or at least resolved greatly by non- 
focusing instruments such as Sw fflQin the near future [^. The fast, accurate localization of 
bursts by non-focusing optics is required to guarantee effective multi-frequency follow- 
up observations, which may be important for understanding the physics of the bursts, as 
demonstrated in recent BeppoSAJ^ observations 

There are many techniques for reconstructing images in non-focusing telescopes. For 
example, in coded-aperture systems, two types of methods are largely used - correla- 
tion and inversion methods. Both methods are subject to either coding noise or quantum 
noise, which leads to spurious fluctuations in images. In order to improve the quality 
of images, a lot of different techniques are also employed such as pixon or maximum 
entropy method (MEM) [0, These methods are computationally intensive and, quite 
often, requires extensive control over somewhat uncertain parameters in order to pro- 
duce satisfactory results. 

More importantly, in the case of applying non-linear refinement algorithms like pixon 
or MEM on images reconstructed by correlation or inversion methods, time-independent 
point spread functions (maybe position dependent) are usually used. Therefore, without 
prior knowledge of the sources intensity change, it is rather difficult to deal with vari- 
able sources when the sources in the FOV have apparent movements due to motions of 
telescope relative to the sky. 

Here we introduce a new method - photon folding - for imaging in non-focusing op- 
tics. The basic idea of photon folding is motivated by the well-known epoch-folding 
method in timing analysis, which is a powerful, yet very simple procedure for finding 
periodicities in data [^^jl^]. The photon-folding method provides a natural way to assess 
the statistical significances of images, which is inherited from the epoch-folding method. 
Our technique uses only projection matrices and can be applied to various types of non- 
focusmg optics. The technique is flexible enough to include corrections for aspect errors 
and non-uniform background. 

^The Constellation-X Observatory is a next-generation X-ray telescope satellite, planned as four individual 
X-ray space telescopes operating together. The Hard X-ray Telescope (HXT) on Constellation-X is designed 
to image X-rays from 5 to 100 keV with resolution better than 1' over a 8' field of view. X-rays are collected 
by graded multilayer reflective optics and imaged with a position sensitive X-ray detector. 

*The Swift is a gamma-ray burst explorer, which will carry two X-ray telescopes and one UV/optical 
telescope to enable the most detailed observations of gamma-ray bursts to date. 

^BeppoSAX is an Italian X-ray astronomy satellite with the wide spectral coverage ranging from 0.1 to 200 
keV. 
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The criterion given by the photon folding technique allows a recursive method for im- 
proving images without human supervision of any parameters, and thus it can be imple- 
mented as on-line data processing. In practice, the relative orientation between telescopes 
and astrophysical sources continues to change during observations and the intensity of 
sources may also vary. In such cases, recursive methods cannot always effectively extract 
the location or intensity of a certain source particularly in the presence of highly variable 
stronger sources. High-order photon folding provides a systematic correction for coding 
noises in such cases beyond what conventional methods can provide. 

We demonstrate the idea using a two-dimensional coded-aperture system. The tech- 
nique can be readily generalized to other systems such as modulation collimators. 



2 Photon Folding Theory 

2.1 Basic Concept 

Consider a typical coded-aperture telescope with a uniform background, where projec- 
tion matrix Pij can have either (closed) or 1 (open mask element). A source Sj from 
direction ; in the sky can generate counts hi at bin i in the detector as 

hi = Y^PijSj + bi, (1) 

/ 

where bi represents the background count at bin i in the detector (H = Lhi). We use 
index S5rmbol / for detector bins and all the other index symbols (;,fc,/) represent sky 
coordinates. Usually the expected image Sj is reconstructed as 

Sj = Y^Gjihi, (2) 

i 

where Gji is a reconstruction matrix, usually given by inversion or correlation methods. 

The epoch-folding method in timing analysis searches for any deviation in the data 
from a flat distribution. Like the x^-test in the epoch-folding method, we can define as 

e = L^, (3) 

where h is the average of hi. In a typical coded-aperture system whose overall trans- 
parency is close to 0.5, from a point source has the largest value while from a flat 
field observation is close to a minimum. Consequently, the above indicates the pres- 
ence of sources within the field of view (FOV), but it cannot reveal the location of the 
sources. 
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Figure 1: Illustration of photon folding: for a given direction in the sky, the detected 
photons are folded into two regions - shadowed/ unshadowed. 



In general, the fraction (p) of open mask elements seen by the detector depends on 
the direction in the sky, i.e. 

where M is the number of detector pixels. We can also define cty, the fraction of detected 
photons which can originate from direction j in the sky, 

/:P/y=l I 

For a given direction in the sky at any moment, one can divide the detector space into 
two separate regions; the region that can have source photons from the given direction 
and the other region that is shadowed by the mask (Fig. ||). If there is no source within 
the FOV, the ratio of the total counts in these two regions will be simply the ratio of 
open/closed mask elements, i.e. aj ~ pj. If there is an appreciable source located at the 
given direction, ocj > pj. The difference of aj and pj indicates the probability of having a 
source at the given direction in the sky. 

For the direction ; in the sky, by folding detected photons into the two regions (Fig. |T]), 
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we define £,j of the photon folding like the x^-test in the epoch folding method as 



^2 _ {Haj-HpjY ^ {H{l-aj)-Hil-pj)y 



Hpj H(l - Pj) 

{(X; — pi)^ 

= w)t^ — (6) 

(1-P;)P; 

In the absence of sources, if Hctj and H(l — ctj) are governed by a normal distribution, 
the of the photon folding follows X2-1 statistics. When pi is the probability density of 

statistics with one degree of freedom, the probability (Qi) of having greater than r\ 
is given as 

Qi(x?>n)=/ vi{)^)d^, (7) 

Jt] 

and then, the confidence level C is related to Qi as, 

l-^ = KQi{x\>Vi), (8) 

where K is the number of sky bins. 

If there is only one source {Sj) within the FOV, the reconstructed image (sy) can be 
directly calculated from otj as 

M{l-pj)pj-''' 

The above formula for sy by photon folding is similar to unbiased balanced correlation 
[11] and becomes identical to the cross correlation method [12] when P,y can have only 



or 1 and there is only one time bin. For the uniformly redundant array (URA) systems 
where pj = p, the above formula reduces to the balanced correlation method where Gji 
is given by 

Gn = frf- (10) 

This is not surprising since photon folding searches for a deviation from a flat back- 
ground distribution similar to the epoch folding method, and the balanced correlation is 
optimized to cancel out a flat background in the detector. 

The immediate advantage of the photon folding is that photon folding provides the 
statistical significance of signals in images without additional calculations of variances. 
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That is, from Eq. @ and (g). 



r2 '■I 
5 jl 



(11) 



M2(l - pj)pj H' 



which is similar to the case of focusing optics where the strength of signals in images 
directly indicates the significance of the signals. Now we derive a generalized photon 
folding formula which can account for aspect errors and non-uniform background. 

2.2 Generalization 

For a given direction (/) in the sky and for a given time bin (t), the detector space can 
be ranked by the effective transparency through collimators, mask elements, detector 
efficiency, etc. One can assign this transparency to the projection matrix, which can now 
have any value from to 1. Thus, in general, the projection matrix Pjy can be defined as 



where 5 is a boundary threshold between shadowed/ unshadowed regions. Let K, M and 
T be the total number of bins in sky, detector and time bin respectively (T could be the 
total duration of the observation when the duration of time bins varies). 

Now, we assume that the detector background is known to follow a pattern fo,(f), 
which can be estimated by flat-field observations or on-source observation of weak sources 
The detected counts are given as 



where bi{t) is the actual background counts at detector bin i and at time bin f. 

For an arbitrary quantity Xi{t) in time bin t and detector bin i, we define the following 
quantities. 



Pij{t) =0...5...1, 



(12) 



hi{t)=l^Pij{t)sj{t)+bi{t), 



H{t)=l^hi{t), 



(13) 



X{t) = '£xi{t), 




1 



1 



(14) 
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where the prime (') notation is used to represent a complementary region. For example, 

Ajt.l[Xl{t),t]+Afk:l[Xi{t),t] = Ak:l[Xi{t),t]. (15) 

We introduce A to simplify the following definitions of a, (5, and p and their explicit 
definitions are given in the appendix. In terms of A, we define 



The above notation allows an intuitive description of the complicated calculations. It 
should be noted that here ^...(0 and p...(0 can be pre-programmed or pre-calculated. 
Again in order to simplify the notation, we introduce the sum on time bins as 

{x} = £x(0. (17) 

t 

Then, the £,j of photon folding is defined as 



_ {H}{H{cKj-^j)}^ 
{H^j}{Hil-Pj)}- 

The first-order reconstructed image sy is given as 

{Sj} _ {H{aj-I3j)] 



{Mipj.,j-p.,jl3j)} 



(18) 



(19) 



In the above formula the aspect errors are handled by using a proper Pij{t) at each time 
bin and the non-uniformity is handled by |3y(f). The formula does not require an estima- 
tion of the overall background level to remove the background pattern. One can express 
the reconstructed image Sj in terms of the true signal counts Sj as 

{fy} ^ {SjiPj-.j-p-.jl^j)} ^ y {Sk{pj:k-P:kl5j)} ^^Q) 

T {Pj:j-P:jl3j] 4- {Pri- P--i^i] 
The second term in the above formula is the coding noise of the system. 
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3 Refinements 



Here we study two refinement procedures of photon folding for suppressing coding 
noise. The first technique, recursive folding, is similar to conventional recursive meth- 
ods (e.g. IROS) [0], which successively remove the side lobes - coding noise of strong 
sources. The second technique, high-order folding, represents a novel approach to the 
coding-noise problem. 

3.1 Recursive Folding 

A simple way to remove coding noise is to photon fold recursively. For a given confidence 
C, the £,j of photon folding will determine detections of signals in the image based on 

Eq. Since the signal Sy is expected to generate counts hj, we replace hi by 

hi <— hi — hi, where 

= j}Prj}/Z (21) 

and repeat the procedure until no excess £,j is found for the given confidence. The final 
image is the sum of Sj at all the recursion steps and the residual. This is very similar 
to many other recursive methods. Recursive folding is deterministic since it does not 
require an initial guess and has a clear stopping point given by the confidence level. This 
feature is missing in some image refinement procedures such as MEM. 

The above formula assumes that the source intensity and detector background level 
are constant during the observation, or that the detector orientation relative to the sky is 
fixed. In practice, this is rarely the case. Without information on source intensity /background 
history, the recursive method does not provide an effective correction to coding noise. In 
this sense, the recursive technique is not a true solution for coding noise. Now we study 
another approach to the coding noise problem. 

3.2 High-Order Folding 

Consider the case where there is only one strong source (s/) within the FOV. Given the 
detection of this source by regular photon folding, one can again apply photon folding 
only in the detector space shadowed by the mask from direction I in the sky. Since in the 
shadowed region there is no contribution from source photons at direction /, this second- 
order photon folding will provide an image which is free of the coding noise from the 
source in direction I. 
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We redefine the total counts and the total number of detector bins only in the shad- 
owed region for the given sky direction / as following: 



H,,{t) = H{t)ar{t), 

MKO= I 1- (22) 



i:Pil{t)<6 



£,f,j for the second-order photon folding will be given by 



2 2 
~ {Hyfiyj} [Hy^yj,] 

_ {Hy]{Hy{0Cyj-^yj)Y 
{Hy^yj]{Hy{l-^yj)y 

The image by the second-order photon folding will be given as 

{S/';};^/ ^ {Hy{0Cyj- ^yj)] 
T {My{pyj.,j- py.,jj3yj)}' 



(23) 



(24) 



The above two formulae are exactly the same as Eq. ( [T^ ) and ( |T9| ) except for an additional 
index I' . The complete expression for an arbitrary order of photon folding can be written 
by simply adding the additional indices for each term (refer to Appendix). In terms of 
true signal counts Sj, 

{Sl'jjjj^l _ {Sj{Pl'j:j - Pl'-.jl^l'j)} _^ y {Sk{pl'i:k - Pl'-.kl^l'j)} 
T {pVy.j- PV:jP>Vj] kj^jj {Pl'j-.j- Pl'-.jlSl'j} 

{Sl{pl'j:l- Pl':ll3yj)} 

H f ^ ^ • (25) 

{Pl'j:j - PV:jPVj) 

The second and third term in the above formula represent the coding noise of the 
system from second-order photon folding. In general, p/'. in the third term keeps the 
coding noise from the strong source S/ very small. For example, in a system with Pij being 
either or 1, there is no coding noise from S/ since py,,,;i = 0. When there is a strong point 
source within the FOV, the above method will provide the best image of weak sources 
within the FOV regardless of the time dependence of the intensity of the strong source 
relative to the background level. 

There are a few problems in applying the above formula in practical situations. First, 
although the second-order folding removes the coding noise from the primary strong 
sources, the coding noise from other weak sources is larger in the second-order folding 
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than in the regular photon folding since the second-order folding utilizes only part of the 
mask pattern or detector space. Such increase of coding noise from weak sources does 
not guarantee the reduction of the overall coding noises in high-order photon folding. 

Second, the above formula might not be suitable for multiple strong sources or an ex- 
tended object. In the presence of multiple strong sources, one might have to rely on third 
or higher-order folding involving the shadowed region for all the directions of strong 
sources. In typical coded aperture systems, p ~ 0.5, so that there is only p" fraction 
of the detector area for the completely shadowed region from n point sources. That is, 
the shadowed region for high-order folding runs out quickly as the number of sources 
increases. 

In order to overcome larger coding noise from weak sources in second-order photon 
folding, the other region - the unshadowed region {{sjj}) by the primary strong sources 
- should be utilized as well as the shadowed region ({s//y}). If the detector background is 
uniform, the second-order photon folding in the unshadowed region will provide similar 
results as in the shadowed region. For example, in the case of the uniform background 
system with being either or 1, the term in {s;y}y^;, equivalent to the third term in 
Eq. (H), is proportional to p;y.; — PulSij which vanishes. 

In general, the images by second-order folding, {s/^jy-^; and {sj'jjj^i will have sub- 
stantially reduced contributions from the coding noise of primary strong sources at di- 
rection I. Both images, however, have larger coding noises from the secondary sources. 
The best second-order photon-folding image will be a linear combination of these two 
images, i.e. 

Y{svj} + {1-Y){sij}. (26) 

The optimal y strongly depends on 5, Pij, and etc. For a given detection of a strong 
source by the regular photon folding, the optimal y for the second-order photon folding 
can be estimated by using only the strong source with constant intensity. In other words, 
we first calculate {s|',y},{s|'y} from hi instead of hi, where 

hi = {{si}Pi!}/T, (27) 

and then use the y in Eq. (^6|), which minimizes 

y{4.} + {l-j){s^.}. (28) 

This procedure is valid when S; is much stronger than Sj, which is the assumption of this 
technique. Here once again we assume the constancy of the source intensity in Eq. (^), 
but due to the terms like pi',„:i in the coding noise from strong source S/, the changes of the 
source intensity are less serious in the second-order folding than in the recursive method. 
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Now, we consider the case of multiple strong sources. When {s;} ^ {sj}, the recon- 
structed images by the regular and second-order folding can be summarized as 

{Sj} = 0{rj){sj} + l^0{ne){sk}, 
y{svi\ + (1 - y){si]\ = 0{rj){sj} + £ 0(T,e){s,} + 0{r,e^){si}, (29) 

where t represents the temporal uniformity of the source with respect to the background 
level and e is the size of coding noise. Since the contributions of source sj. (the second 
term in the right hand side of the above formula) should be similar in both images from 
the first-order and second-order folding, we expect that the following formula provides 
a more reliable estimate for Sj in the case of multiple strong sources. 

£ (y{s,,^} + {l-r){Sij})-{L-l){Sj} 

= 0{rj){sj} + Y^O{ne^){s,}, (30) 

where L is the number of sources detected in the first-order photon folding and y is 
calculated by minimizing 

X r{s?.,} + (i-r)H}. (31) 

The left hand side of Eq. ( |30| ) is expected to reduce the coding noise below that from 
the first-order image and is less sensitive to the intensity variation of the strong sources 
compared to the recursive method. 

4 Simulations 

We demonstrate the photon folding technique on a coded-aperture imaging system. The 
system consists of 64 x 64 random mask elements {Pij= or 1) and the detector has 32 x 
32 elements. 

The first example has four point sources and one extended source within the FOV. 
Two of the point sources are very strong and the other sources are relatively weak. During 
the simulated observation we assume that the relative orientation of the mask pattern 
with respect to the sky has rotated 90 degrees from the first half to the second half of the 
measurement. This is a somewhat extreme case of the change of telescope orientation 
relative to the sky during the observation. When reconstructing images, we combine all 
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the data and we assume that the intensity change of the sources during the observation 
is unknown, i.e. the intensity is assumed to be constant. The aspect system and spatial 
resolution of the detector are assumed to be perfect. 

The upper plots in Fig. ^ show the true sky image and the simulated detector counts. 
We assume that the detector background exhibits a quadratic spatial dependence, result- 
ing in an enhancement of counts at the edges. Such a pattern is somewhat common and 
we assume that the pattern is known by previous flat field observations or calibrations. 

In Fig. ^ the intensity of sources in the sky did not vary during the simulated obser- 
vation. The lower three plots show the reconstructed images by regular photon folding, 
recursive folding, and second-order folding. The regular photon folding, which is equiv- 
alent to the cross correlation method in this case, successfully detects the strong sources, 
but the coding noise of the strong sources buries features of weak sources. The recursive 
folding presents the best image in this case, and its quality is limited only by random 
Poisson noise. Second-order folding also reconstructs a decent image which reveals the 
fine structure of weak sources. The noise in the image by second-order folding is mostly 
random Poisson noise and coding noise of weak sources. 

Now, in Fig. ^ we have an extreme case of a more realistic situation. The intensity of 
two strong sources changed dramatically from the first to the second half of the obser- 
vation. If the relative orientation between the telescope and the sky did not change, the 
results would be similar to those in Fig. ^ But here we assume that the relative orienta- 
tion changed 90 degrees. 

Under the assumption of being imaware of the source intensity changes, the bottom 
plots in Fig. |^ show the three reconstructed images. The recursive method fails to show 
the fine structure of weak sources. It should be noted that this result is generally true 
for any type of recursive method without prior information on source intensity changes. 
On the other hand, the second-order folding produces almost the same image as in the 
previous case, i.e. only limited by the random Poisson noise and coding noise of weak 
sources 

In order to see the general performance of the recursive method and the second-order 
photon folding, we calculate SNRs of a few simple simulations with a flat background 
pattern. In the following figure, we show the SNR of a point source in four distinct 
situations (1 cr distribution of the simulation results). Each observation consists of two 
measurements as in the previous example (90 degree offset of the relative orientation). 

Fig. ^ (a) shows the case of a steady single source within the FOV. The second-order 
folding and recursive folding generate images with the maximally allowed SNR, while 
SNRs from regular folding are limited by coding-noise. In Fig. ^ (b), the source intensity 
dropped to zero for the second half. While second-order folding performs similarly to 
the previous cases, the performance of recursive folding is limited by coding noises. 

In Fig. ^ (c), there are five steady point sources within the FOV. In this case, the re- 
cursive folding produces perfect images, while second-order folding shows its limitation 
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Figure 2: Simulation of steady sources: the upper plots are the simulated sky image, total 
counts, and background counts in the detector. The lower plots show the reconstructed 
images by regular photon folding, recursive folding, and second-order folding. The ori- 
entation of the telescope relative to the sky has changed 90 degrees from the first to the 
second half of the observation. 



due to incomplete correction of the coding noise from multiple sources. In Fig. ^ (d), 
there are five variable point sources, and the intensity of three sources drops from max- 
imum to zero and the intensity of the other two rises from zero to maximum from the 
first to the second half of the observation. It is clear that if there is a change in source 
intensity, second-order folding is the better choice, while in the case of many multiple 
steady sources recursive folding is the optimal method. 



5 Discussion 



The photon folding and associated techniques can be applied in many non-focusing in- 
struments. Recursive folding and second-order folding are complementary to each other 
in various realistic situations. Here we discuss some of the fundamental issues. 
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5.1 Computational Issues 



Computational time for regular photon folding is similar to that of the cross correlation 
method. Since the photon folding automatically takes care of aspect errors and non- 
uniform background, the overall computational burden for photon folding is similar to 
regular correlation or inversion techniques. 

If we let N be the number of calculations for photon folding, the recursive folding re- 
quires NrN calculations where Nr is the total number of recursion steps. depends on 
intensity distribution of sources in the FOV. Second-order folding requires {2Ns + 1)N 
calculations, where Ng is the number of strong sources detected in the regular photon 
folding. The extra N comes from the estimation of the optimum y for the given detec- 
tion. Since photon folding does not involve inversion of matrices, this number of cal- 



IstHalf of 
True Sky Image 



2nd Half of 
True Sky Image 




Figure 3: Simulation of variable sources: the upper plots show the true sky image in the 
first and second half of the observation. The orientation of the telescope relative to the 
sky has changed 90 degrees from the first to the second half. Reconstructed images by the 
regular photon folding, recursive folding and second-order photon folding are shown in 
the lower plots. 
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Figure 4: SNRs in various cases: within the FOV, there is (a) a steady point source (b) a 
variable point source, (c) five steady point sources, and (d) five variable point sources. 
The overall intensity of each source in multiple-source cases are the same. 



culations is not a problem even for a huge number of sky or detector pixels due to the 
improvements of modem computing power. 

Both recursive and second-order folding can be implemented on-line. While the re- 
cursive foldiag is somewhat straightforward in on-line processing, second-order folding 
requires more caution. The successful image reconstruction of second-order folding de- 
pends on limiting the number (L) of sky pixels for the first-order detection and yet includ- 
ing enough of them to remove most of the coding noise. Selecting a few of the strongest 
sources for the first detection would be adequate for an automatic implementation. 
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Figure 5: Recursive vs Second-order photon folding: (a) a steady source (b) a variable 
source, (c) five steady sources, and (d) five variable sources. The x-axis is the SNR by 
second-order photon folding relative to the simulated SNR, and y-axis the SNR by recur- 
sive folding folding relative to the simulated SNR. Each shaded region represents each of 
the four cases of simulations in Fig. § where simulated SNR is between 50 and 1000. 



5.2 Recursive vs Second-order folding 

These two techniques are somewhat complementary to each other in various situations. 
Fig. ^ summarizes the simulations in the previous section (Fig. ^) in terms of relative 
SNR. In Fig. ^ each of the shaded regions represents each case of the simulations in 
Fig. § where simulated SNR is between 50 and 1000. For a given value of the simulated 
SNR, recursive folding produces better results in the case of many steady sources, while 
second-order folding does so in the case of variable sources. 

The second-order photon folding is less robust against aspect errors or non-uniformity 
of the background than the recursive folding. Having uncorrectable aspect errors is 
somewhat equivalent to having multiple sources in the FOV. The non-uniform back- 
ground in the second-order image ({s/y}) from the unshadowed region cannot cancel 
out in photon folding, depending on the strength of the strong source S/ relative to the 
background level. 

5.3 Other variations of photon folding 

So far we use only two folded bins - shadowed/ unshadowed region - for photon fold- 
ing. In general, increasing the number of folded bins does not boost SNR. However it 
may be useful to have more than two folded bins when effective transparencies have 
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more than two distinct values or there are severe non-uniformities in the detector back- 
ground. One can also utilize higher-order folding beyond second-order, such as S/')t';/ 
Siikj, etc. Although higher-order foldings are computationally intensive, the use of higher- 
order folding could remove the limitation of the second-order folding method in multiple 
strong source cases. Further studies are required to find a proper form of higher-order 
folding to suppress even more of the coding noise and also new types of optics scheme 
can be designed to optimize for high-order photon-folding image reconstruction. 



5.4 Applications 

The photon folding technique is very versatile, and it can be applied to many different 
types of experiments. Applying the photon-folding method to other non-focusing sys- 
tems like Fourier grid systems is quite straightforward. In the case of modulation colli- 
mator systems, the location of the source is identified by the temporal modulation rather 
than the spatial modulation in the detector (usually detectors for modulation collimator 
systems do not require spatial resolution). 

For example, consider a typical rotation modulation collimator. For a given direction 
in the sky, the fraction (p) of the shadowed area in the detector changes with time. For 
a given unit rotation and a given direction in the sky, one can rearrange count rates in 
the order of the fraction (p) instead of time, and then apply a folding procedure with 
a necessary amount of birming. This is somewhat similar to regular epoch-folding in 
timing analysis. The signal from a point source usually fluctuates between a maximum 
and minimum during the unit rotation. The temporal resolution, relatively finer than 
the spatial resolution of a typical detector, might allow effective usage of more than two 
folded-bins in photon folding. 

In coded-aperture systems, photon folding may allow use of non-URA mask patterns. 
For example, EXIST is a coded-aperture, wide field of view survey mission, with a wide 
energy range (~ 10 - 600 keV) [^. Curved mask patterns are being studied for EXIST 
to utilize maximal available FOV without substantial collimators between the mask and 
detectors. In order to cover the wide energy range (10 - 600 keV) in EXIST, two scale 



mask patterns can be used with energy dependent transparency [[13[]. Such mask patterns 
may provide additional advantages over conventional mask patterns, but the optimal 
imaging reconstruction scheme is not yet available. The photon folding technique may 
be useful for non-conventional mask patterns. Particularly second-order folding is very 
interesting for EXIST. The high sensitivity of EXIST is likely to result in detection of 
multiple, possibly variable sources in each field, and second-order folding would allow 
imaging of the weak sources. 
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6 Conclusion 

A new imaging technique - photon folding - is introduced for non-focusing telescopes. 
The technique is quite robust against common imaging problems like aspect errors and 
non-uniform background. Its performance is demonstrated by a two-dimensional coded- 
aperture system and photon folding can be applied to other types of imaging telescopes. 
Two refinements of photon folding are presented - recursive and second-order folding. 
In particular second-order photon folding is suitable for imaging weak sources in the 
presence of highly variable strong sources regardless of the changes in telescope orienta- 
tion relative to the sky. 
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Appendix A 

For the regular photon folding. 




(32) 



which reduces to Eq. when Pij is either or 1 (constant). 




(33) 



where = T.bi{t). For a flat background pattern, ^y(t) = Pj;{t) 
And 




(34) 



If Pij{t) is either or 1, p.y(f) = pj.j{t) = pj.{t). 
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Now, for second-order photon folding, 

A}j:[hi{t),t] 



Likewise, 



where 



aij(t) = 



Fij{t)>S Pij(t)>6 



•^IV) i:Pu{t)>S, 
Pij(t)>S 

Pi/(t)>* 

i:Piiit)>5 



Biit) = B{t)l3i{t), 
Mi{t)= £ 1. 

i:Pn{t)>5 



(35) 



(36) 



(37) 



In ociij{t), ^iij{t), Pi'j:{t), Pi'j;j{t) and pi':j{t), the summation condition, Pii{t) > 5, in the 
above equations will be replaced by Pii{t) < 5. Therefore, they satisfy the foUowings 
relations. 



I5lj{t) + Pijit) = 1. 



(38) 



The complete expression for an arbitrary order of photon folding can be written sim- 
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ply by keeping additional indices in front of each term. 



i^h-in {H,^.../„/3;,.../„/}{H;^...;„(l-^/,..w)}' 



(39) 
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